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g : 

In present paper we suggest a new universal approach to study complex systems by microscopic, 

mesoscopic and macroscopic methods. We discuss new possibilities of extracting information on non- 

stationarity, unsteadiness and non-Markovity of discrete stochastic processes in complex systems. 

We consider statistical properties of the fast, intermediate and slow components of the investigated 

processes in complex systems within the framework of microscopic, mesoscopic and macroscopic 

CZ3 \ approaches separately. Among them theoretical analysis is carried out by means of local noisy 

* i-H 

time-dependent parameters and the conception of a quasi-Brownian particle (QBP) (mesoscopic 

-i— j , 

approach) as well as the use of wavelet transformation of the initial row time series. As a con- 
crete example we examine the seismic time series data for strong and weak earthquakes in Turkey 
(1998, 1999) in detail, as well as technogenic explosions. We propose a new way of possible solution 
Q ■ to the problem of forecasting strong earthquakes forecasting. Besides we have found out that an 

i ^ i ■ unexpected restoration of the first two local noisy parameters in weak earthquakes and technogenic 

I ' explosions is determined by exponential law. In this paper we have also carried out the compar- 

ison and have discussed the received time dependence of the local parameters for various seismic 
(f) ' phenomena. 

\o ■ 
o ■ 

1 I. INTRODUCTION 

o : 

Nonstationarity, unsteadiness and non-Markovity are the most common essential peculiarities of stochastic pro- 
^ , cesses in nature. The existence of the similar properties creates significant difficulties for the theoretical analysis of 
real complex systems Q. At present, methods, connected with localization of registered or calculated parameters 
for the quantitative account of the dramatic changes caused by the fast alternation of the behavior modes and in- 
^ \ termittency, came into use. For example, the time behavior of the local (scale) Hurst exponents was defined in the 
recent work of Stanley et al. to study multifractal cascades in heartbeat dynamics [lj and to analyze and forecast 
earthquakes and technogenic explosions in Ref. 2]. The application of the local characteristics allows one to avoid 

_ _ i 

difficulties connected with nonergodicity of the investigated system and gives a possibility to extract additional valu- 
able information on the hidden properties of real complex systems. From the physical point of view, this approach 
resembles the use of nonlinear equations of generalized hydrodynamics with the local time behavior of hydrodynamical 
and thermodynamical parameters and characteristics. 

It is well known that one of the major problems of seismology is to predict the beginning of the main shock. Though 
science seems still to be far from the guaranteed decision of this problem there exist some interesting approaches based 
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on the peculiar properties of precursory phenomena [il 0, 0, Q, S, Q ■ Another important problem is recognition and 
differentiation of weak earthquakes and technogenic underground explosions signals. One of useful means of solving 
this problem is the defining of corresponding local characteristics Q . 

In the present work we suggest a new universal description of real complex systems by means of the microscopic, 
mesoscopic and macroscopic methods. We start with a macroscopic approach based on the kinetic theory of discrete 



stochastic processes and the hierarchy of the chain of finite-difference kinetic equations for the discrete time correlation 
function (TCF) and memory functions 

The mesoscopic phenomena of the so-called "soft matter" physics, embracing a diverse range of system including 
liquid crystals, colloids, and biomembranes, generally involve some form of coupling of different characteristic timc- 
and length-scales. Computational modelling of such multi-scale effects requires a new methodology applicable beyond 
the realm of traditional techniques such as ab initio and classical molecular dynamics (the methods of choice in the 
microscopic regime), as well as phase field modelling or the lattice-Boltzmann method (usually concerned with the 
macroscopic regime). As for complex systems, we propose to consider intermediate and slow processes within a unified 
framework of mesoscopic approach: by means of local time behavior of the local relaxation and kinetic parameters, 
local non-Markovity parameters and so on. For this purpose we introduce the notion of quasi-Brownian motion in a 
complex system by coarse-grained averaging of the initial time series on the basis of wavelet transformation. 

As an example we consider here the local properties of relaxation or noise parameters for the analysis of seismic 
phenomena such as earthquakes and technogenic explosions. The layout of the paper is as follows. In Sec. ITT1 we 
describe in brief the stochastic dynamics of time correlation in complex systems by means of discrete non-Markov 
kinetic equations. Basic equations used for these calculations are presented here. The local noise parameters arc 
defined in Sec. IIIII Section llVl contains results obtained by the local noise parameter procedure for the case of seismic 
signals. The models of the time dependence of the local parameters are given in Sec. [V] The basic conclusions are 
discussed in Sec. IVII 



II. THE BASIC DEFINITIONS IN KINETIC DESCRIPTION OF DISCRETE STOCHASTIC 

PROCESSES 

A. MACROSCOPIC DESCRIPTION IN THE ANALYSIS OF STOCHASTIC PROCESSES 

A lot of different existent processes, such as economical, metheorogical, gravimetrical and other, are registered as 
discrete random series Xi of some variable X. This random variable X can be written as an array of its values 

X = {x(T), x(T + t), x{T + 2t), x(T + kr), x(T + (N - l)r)}. (1) 

Here a time step (or a time interval) r is a constant, T is the time when the registration of the signal begins, (N — l)r 
is the duration of the signal detection. 



The average value (x) and fluctuations Sxj are denned by the following expressions, correspondingly 

1 N ^ 

W = ]vE a: ( T + J T )' 6x j =5x(T + jT)=x(T + j T )-{x}. (2) 

3=0 

From the fluctuations of the considered random variable Sxj we can form fc-component state vector of the system's 
correlation state 

A° = A^(O) = (5x(T), Sx(T + t), 5x(T + (k - l)r)) = (Sx , 5x u Sx k -i). (3) 

The time dependence of the correlation state vector A can be represented as a discrete m-step shift 

A™ +fc = A£+ fc (t) = (5x(T + mT), Sx(T + (m + l)r), <5x(T + (m + fc + l)r)) = (fe , <fci, fom+fc-i). (4) 

Then by analogy with the papers we can write the following normalized TCF 

= (A° f _ 1 _ m A^_ 1 ) = (A° f _ 1 _ m (0)A^_ 1 (t)) = (4-i-mg = T + mr, r)A° f _ 1 _ m (0)) 

(Ajv-l-mAjv-l— m) |AAr_i_ m (0)| 2 I A.jy_ 1 _ m (0)| 

where angular brackets indicate the scalar product of the two state vectors. On the other hand, the time dependence 
of the vector A r ^_ 1 _ m (T + t), t = mr, can be represented formally with the help of the evolution operator U (t', t) as 
follows: 

A^_ X (T + t) = U(T + mr, T)A° r _ 1 _ m (T) - U(t, 0)A° N _ 1 _ m (0). (6) 

The last one has the property: U{t,t) = 1. Actually, one can write down formal discrete equation of motion with the 
use of the evolution operator U (t',t ) (see Appendix A for more details). 

It was shown in Refs. [A IllL Il2| that the finite-difference kinetic equation of a non-Marko v ty pe for TCF Mo(t) 
can be written by means of the technique of projection operators of Zwanzig'-Mori's type [lolllq as 

AM (t) ra_1 

A ° V ' = AiMo(i) - t-Ai J2 Mi(jr)M (f - jr). (7) 

j=o 

Here the first order memory function M\(jr) appears, Ai is the eigenvalue of Liouville's quasi-operator C and Ai is 
the relaxation noise parameter, which are characteristics of the investigated process. Possible methods of defining 
quasioperator £ are presented in Appendix A [see Eqs. (|A3|) . i|A5l) and (|A7() ]. It should be recorded that Eq. (JJJ is 
the first kinetic finite-difference equation for initial TCF Mo(t). With the use of the same procedure of projection 
operator we can obtain the chain of kinetic finite-difference equations of the following form: 

AM (t) m_1 

g-i^ = A i M i _ 1 (t)-rA i ^M i (ir)M i _ 1 (t-ir), i= 1,2,3,.... (8) 

3=0 

Here Mi(Jr) is the memory function of the ith order, whereas A^ and Ai are noise parameters: 



, _ . (Wn-tlWn-i) _ . (W n ^CW n ) 



Here W„ are the dynamical orthogonal variables, obtained by the Gram-Schmidt orthogonalization procedure 

(W„,W m ) = 5 n , m (|W r , 



•«l 2 }, 
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where 8 n>m is the Kronecker's symbol, 

W =A°(0), Wi = [i^-Ai]W , 



Wa = [t£-Aa]Wi-AiWo,.... (10) 

From Eq. I|10(l it is obvious that in the cited Gram-Schmidt procedure from each new vector of state one should 
subtract the projection on all the previous vectors. Thereafter the orthogonalization (|10|l is complete. 

A chain of integro-differential equations JSJl arise as a result of the use of projection operator technique to define 
different correlation functions in physical problems (for example, TCF of density fluctuation in Inelastic 

Neutron Scattering HQ and Light Scattering investigations, velocity autocorrelation function and others can be 
received in the specified way). However, the quest for the physically based way of closing the chain of equations 
and finding the spectra of the initial TCF are essential moments in these challenges. Here the situation is different. 
Namely, the initial TCF can be calculated directly from the experimental data. Then memory functions Mi(t) and 
noise parameters A.;, A, are similarly calculated from the experiment. All these functions and parameters make it 
possible to carry out the detailed analysis of the random process. The noise parameters Xi and Aj are the relaxation 
characteristics of the experimental time series, which contain information of various modes passing and changing. The 
macroscopic approach presented above is based on the calculation of memory functions, power spectra, dynamic vari- 
ables, relaxation parameters. It suggests the investigation of the system as a single whole. The global characteristics 
calculated on the basis of all the time series contain hidden information about various modes of the system behavior. 
As a rule, this information is difficult to extract and analyze. For this reason it is necessary to develop a mesoscopic 
description and introduce a local noise and relaxation parameters Xi and A^. 

B. MESOSCOPIC DESCRIPTION IN ANALYSIS OF STOCHASTIC PROCESSES 

It is well known, the mesoscopic conception is one of the possible ways to deal with random processes in complex 
systems. It consists of extension of the domain of the dynamical equations mesoscopic variables and of introduction 
of a some local time interval. Such quantity as X (it can be particle position, mass density and so on) is defined 
on the mesoscopic space In addition we introduce number M as its extensive quantity. This number M should 
satisfy the condition: 

1 <C M <C N, (11) 

where N is the length of the initial experimental sampling. Let us take a working window of the fixed length M. By 
superposing this window on the initial sampling X, we choose all elements, incoming into it, as a separate sampling 
£o- Further, let us execute one time step r shift of this working window to the right and obtain another local sampling 
of the length M. Executing this procedure (N — M + l) times, one can obtain the same quantity of the local samplings 
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of length M: 

6 = Zo{x(T),x{T + t),x(T + 2t),...,x(T + (M-1)t)}, 
6 = £i{x(T + t),x(T + 2t),x(T + 3t),...,x{T + Mt)}, 

**•) 

itf-M = £n-m{x(T+(N-M-1)t),x(T+(N-M)t),x(T+(N-M+1)t)...,x(T+(N-1)t)}. (12) 

The obtained local samplings £j form the array, which represents the time dynamics of the investigated process, 

f(0 = 6,6, 6v-m}. (13) 

Then, in accordance with the procedure described in the last subsection All . we can define fluctuations by Eq. J5J), 
calculate TCF with the help of Eq. 10, calculate memory functions and parameters Aj and Aj by Eqs. for every 
sampling £j. However, parameters Aj and Aj will be characteristics of concrete jth sampling only. To characterize the 
local properties of the initial time data X, it is convenient to represent their local time dependence in the following 
way: 

Ai(t') = {A,(T + (M - l)r), Xi(T + Mr), \(T + (M + l)r), . . . , A,(T + (N - M)r)}, 
Ai(t') = {Ai(T+(M - l)r),Ai(T + Mr), A«(T + (M + l)r), . . . , A,(T + (N - M)t)}. (14) 

Operating in the similar way, one can execute cross-over from the macroscopic description of the whole system 
to the mesoscopic one. The offered approach is very convenient for the description and analysis of non-stationary 
stochastic processes. It allows to depart from the global macro-characteristics, which carry only averaged minor 
information about the whole investigated process, and to turn to the stochastic description with the use of the local 
characteristics, and, as a result, to execute a more detailed analysis of various dynamic states of the system. 



C. CONCEPTION OF ONE-DIMENSIONAL DYNAMICS OF QUASI-BROWNIAN PARTICLE 

Let us consider the motion of a large Brownian particle in a dense medium composed of light molecules and restrict 
it by a simple one-dimensional case. The coordinate Xi and the velocity Vi are random variables of a Brownian 
particle. The quantity r represents average time between the two successive collisions of liquid molecules, T is the 
initial moment of time, the variable M characterizes a local size (mass) of a Brownian particle. It is obvious, that a 
Brownian particle must have a larger mass in comparison with liquid particles [M S> 1), therefore it is more inert. 
Then it is convenient to define the coordinate of a Brownian particle at moment t' as an average value of sampling 
£i. For example, we obtain from £o: 

M-l 



x{T) + x(T + t) + x(T + 2t) + . . . + x(T + (M - l)r) _ 1 



M M 

3=0 



= TjJ2 x ( T +^y ( i5 ) 



The quantity yo defines the coordinate of "the center of mass" of a Brownian particle at the initial time moment 
t' = 0. By analogy, it can define the position of a Brownian particle at the next time moment t' = r and so on. As a 



result, we obtain a new time discrete series Y(t') as: 

Y (t') = {Vo, 2/1,2/2, ...,2/jv-m}- (16) 

The velocity of a Brownian particle is 

Vi = . (17) 

T 

So, for example, for the initial velocity vo we obtain from Eqs. I|12[i . I|15|l and l|17(l : 

2/1 - 2/o ^(T 1 + Mr) - a:(T) 

«0 - - 7T ■ (1°) 

r Mr 

Obviously, £/ t/ie larger one is M , the smaller one is the velocity of a Brownian particle. So, we obtain a discrete set 
of velocity values for a Brownian particle at equally small time interval r 

V(t') = {v ,vi,v 2 , ■ . . ,v N - M -i}- (19) 
Eqs. (|T^)l and lfT§|l define time dependence of random variables Y(t') and V(t'). 

1. Generalized Langevin equation with discrete time 

In accordance with Eq. we define the correlation state vector, components of which are the fluctuations of the 
particle position, 

B° k = {5y , Syi, Sy 3 ,..., Sy k ^}. (20) 

Then, by analogy to Eq. Q time dependence of the vector B can be considered a discrete m-step time shift. For the 
vector B the following normalized TCF can be written with the help of Eq. (JSJ 

m = {B h~ mB T } - (2D 

\ ti N-l~7n\ 

The last one describes the time correlation of the two different correlation states of the system. 

Now we introduce the linear projection operator in Euclidean space of the state vectors 

B(0))(B(0)B(t)) _ B(0))(B(0) 
QB \BW 2 " (0)> m Q ~ (B(O)B(O))' (22) 

This operator has the necessary property of idempotency Q 2 = Q. The existence of projection operator Q allows to 

introduce the mutually-supplementary projection operator R as follows: 

R=l-Q, R 2 = R, QR = RQ = 0. (23) 

It is necessary to note that the projectors Q and R are both linear and can be recorded for the fulfillment of projection 
operations in the particular Euclidean space. The projection operators Q and R allow one to carry out the splitting 
of Euclidian space of vectors B, where B(0) and B{t) 6 B, into a straight sum of the two mutually supplementary 
subspaces in the following way: 

B = B' + B", B' = QB, B" = RB. (24) 
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Let us consider the finite-difference Liouville's Eq. (|Af I) for the vector of fluctuations of a Brownian particle position 

A B - +fe (t)=i£(t,r)B- +fc (t). (25) 

Affecting the last equation by the operators Q and R successfully, we can split the dynamic equation l|25|l into two 
interconnected equations in the two mutually-supplementary Euclidean subspaces: 

^B'(i) = iQC(Q + R)B{t) = ^£llB , (^) + iC 12 B"(t), (26) 

£iB"(i) = iRC(Q + R)B(t) = il 2l B'{t) + iC 22 B"(t). (27) 
In the above equations the matrix elements £y of the quasi-operator C have been introduced 

£■11 — QCQ, C\ 2 — QCR, L 2 \ = RCQ, C 22 — RCR, 

\jl 2 i L 22 J 

Operators Cij act as follows: Cu - from a subspaces B' to B'; £ 12 - from B" to B'\ £ 21 - from B' to B"; and £22 - 
from B" to B" . 

To simplify the Liouville Eqs. H26|) and (|27|l . we exclude the irrelevant part B"(t) and construct the closed equation 
for the relevant part B'(t). For this purpose, it is necessary to obtain a step-by-step solution of Eq. (|27Jl 

^> = B "" + T) - B " ( " = iLjt® + ^B» (l ), ,29, 

B"(t + r) = B"(t)+zr£ 21 B'(t)+iT£ 22 B"(t) 

= (1 + iT£ 22 )B"(i) + ir£ 2 iB'(t) 

- C/ 22 (i + T,t)B"(t)+ iTifiai^ + T.tJB' (f), (30) 

where U 22 (t + r, t) = 1 + irC 22 (t + r, i) is the operator of a time step shift. 

With the help of Eq. 13U[1 we can derive the following expression for the next time step: 

B"(t + 2r) = U 22 (t + 2T,t + T)B"(t + r) + irC 21 (t + 2r,t + t)B' (t + t) 

= U 22 (t + 2r, t + T)[U 22 (t + r, t)B"(t) + iTC 21 (t + r, t)B'(t)] 

+i-r£ 2 i (f + 2r, t + r)B'(f + r) 

= l7a3(t + 2T,t + T)Cr22(t + T,t)B // (t)+iT[Cr 2 2(t + 2T,t + r) 

x£ 2 i(< + r,t)B'(t) + C 2l {t + 2r,t + r)B'(t + t)]. (31) 

Generalizing this result in case of the mth discrete step we find the following final result 



B"(i + mr) = if J] U 22 (t+(j + l)T,t + jT) \B"(t) 



3=0 

rn—1 I rn — 2 



+iT E W II ^ + C?" + 2)r, t + (f + l)r) 
x£(t + (j + l) r , t + jr)B'(t + jr). (32) 
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Here T denotes the Dyson operator of chronological ordering. Substituting the irrelevant part of Eq. (|26(l in the 
right side of Eq. Ij32|l . we obtain the closed finite-difference equation for the relevant part of the correlation state 
vector 

^B'(t + mr) = iCu(t + (m + l)r,t + mr)B'(t + mr) + itvi{t + (m + l)r,t + mr) 



x 



m — 1 ] m — 1 m—1 

({ f J] U 22 (t + (j + l)r, t + jt) \ B"(t) - r £ {T J] E/ 22 (t + (j' + 2)r, t 

3=0 I 3=0 j>=j 



+(f + l)r)}£ 21 (t + (j + l)r, t + jr)B'(t + jr)). (33) 

Substituting Eqs. I|22[l and 124|) in Eq. I|33|) . we derive a finite-difference kinetic equation of a non-Markov type for 
TCF b(t) 

^ = - rAj ^ Mx(t - jr)6(jr). (34) 

Here the TCF Mi(t) is the first order memory function 



(Wi(0)Wi(t-jr)) 

|W!(0)| 2 



A*i(*-Jt) = ( i\\ i2 . ( 35 ) 



Here Ai is the frequency relaxation parameter of the first order with a square frequency dimension, and is defined by 
Eq. 10. Then if the dynamic variable Wo = B°(0) represents the fluctuations of a Brownian particle position, the 
dynamic variable Wi = i£Wo — AiWp contains fluctuations of pulses of a Brownian particle [see Eqs. I|10|) ]. The 
function M\ (t) is time correlation function of fluctuations of a Brownian particle velocity. 

Defining the corresponding projection operators to new dynamic variable Wi and repeating the above described 
procedure, we find the finite-difference kinetic equation for M\(t) 

= \ 2 M 1 (t) - tA 2 J2 M 2 (t - jrJMx (jt) . (36) 

Here A2 and A2 are the frequency relaxation parameters of the second order, M 2 (i) is the second order memory function 
or memory function of the velocity correlation function for a Brownian particle (memory friction) |l3.l2(i l2lll22| . In 
fact, Eq. (|36[) is a discrete finite-difference generalized Langevin equation (GLE). So, M 2 (t) can be associated with 
TCF of the stochastic Langevin forces, for which the similar equation can be received 

AM (t} m—1 

-^1 = X 3 M 2 (t) - t A 3 A/ 3 (t - jr)M 2 (jt) (37) 

j=o 

with the third order frequency relaxation parameters A3 and A3, and the memory function of the third order M^(t) 
respectively. 

The three Eqs. I)34[l. (|36|l and (|37|l are the exact consequence of microscopic discrete finite-difference equations of 
motion. Calculation of the memory function Mi(t) and the relaxation frequency parameters Xj, A& are the central 
point here. 
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initial s amp ling 

QQHH 



In case of a Brownian particle presented above, dynamic variables Wo and Wi are the position and pulses of 
random particles. In fact, they can be taken as characteristics of any other non-stationary process. The averaging 
operator of sampling with the length M, A = l/M^^J^ 1 [see Eq. <|15|) ]. applied to the intermediate local sampling, 
can be changed by any other operator, depending on the goal of the investigation. The operator A allows one to get 
clear of sharp fluctuations in the initial sampling of data X(i) and to replace it by an other Y(t'): X(t) — > Y(t'), 
which contains the results of coarse-graining averaging. However, it is not always convenient to average a local 
sampling. In these cases the operator A can be replaced by an other operator, which allows one to obtain only one 
number from every local sampling. In general case, an another, a more universal method of transformation of the 
initial s amp lin g int o the sampling with some specified (required) characteristics can become discrete wavelet-transform 
l2Sj . which is defined by the following equation: 

W(j, k)=Y,J2 X{k)2~ j ~ 2 V{2~ j n - k). (38) 
j k 

Here X(k) is the sampling and *&(t) is a time function with fast decay called mother wavelet. The following 
analysis can be applied to the new transformed data W(j, k) according to the algorithm described above. Namely, 
memory functions Mj(i) and frequency relaxation parameters and A* can be calculated for the transformed data. 

2. Analogue of Green-Kubo relation for diffusion coefficient for time discrete system 
According to the theory of random walkers the mean-square displacements of a Brownian particle can be defined 

as 

/+oo 
Ay 2 $ 1 (y,t)dy = 2Dt, (t -> oo). (39) 

-oo 

Here $i(y, i) is the density of probability of being of particle at point y at time moment t. The Eq. (|39J) is a well 
known Einstein relation for continuous displacements. 
According to Eq. I|39|l the displacement during time t is 

Ay* = y(t) - y(0). (40) 

In case of the discrete time Eq. (|40|) can be rewritten in the form 

Ay(t) = Ay(T + jr) = y(T + (j + to)t) - y(T + jr), t = mr, m > 1. (41) 

Then the mean-square displacement of a Brownian particle is 

^ N—m—l 

(Ay 2 (T+(j + m)r)) = - V [y(T + (j + m)r) - y(T + jr)} 2 , N > m > 1, (42) 

N — to 

3=0 

where N — m is the quantity of " possible ways" . The diffusion coefficient takes the following form 

N-m—i 

2m(iV — m)i 



D = - — — V [y(T+(j + TO)r)-t/(T + jr)] 2 , at > to > 1 (or t -> oo ). (43) 

/ —TOT * — ' 



3=0 
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Let's consider separately the sum in the last expression in terms of the velocity [see Eq. <|17|) ] 



N-m-l 



t y (T+(j + m)T) -y{T + j T )} 2 



3=0 



= [y(T + mr) - y(T)] 2 + [y(T + (1 + m)r) - y(T + r)} 2 + . . . + [y(T + (N - l)r) - y(T + (N - m - l)r)] 2 

N „ ' 



(N — m) of square brackets [...] 



= [y(T + mr) - y(T + (m - l)r) + y(T + (m - 1)t) - ... - y(T)} 2 + ... 



[y(T + (N - 1)t) - y(T + (N - 2)r) + tj(T + (JV - 2)r) - ... - y(T + (AT - m - 1)t)] 2 
= r 2 [i;(T + (m - l)r) + «(T + (m - 2)t) + ... + «(T)] 2 + . . . 
+ t 2 [v{T + {N- 2)t) + v(T + (N- 3)r) + . . . + v(T + (N - m - l)r)] 2 



N — m~l m—l+j 




j=0 j=0 



Then the expression for the diffusion coefficient (|43|l can be rewritten as: 



N — m—l m—l+j 



D = 



T 



£ ]T v(T + kr) 



(45) 



2m(N — to) 



J=0 k=j 



Eq. (|45|l is the discrete finite-difference analog of the famous Green-Kubo relation for the diffusion coefficient. 
Given relation has been obtained from Einstein equation l|39JI for a discrete system. The asymptotic limit t — > oo can 
be replaced here by the similar condition N — > oo, m — > oo, N > to. 



Seismic data represent discrete random series, which is recording of displacements of the Earth's surface. There- 
fore, we can use the above-stated formalism to analyze seismic data. In particular, the local dependence of various 
can be serve as additional source of information on properties of objects. The noise parameters 
Aj and A, are very sensitive to the presence of a nonrandom component in the sampling. The change of the character 
of the correlated noise and, the appearance of the additional signal in the sampling can cause the alternation of these 
parameters. So, the time behavior of the local parameters Aj, is important and informative for analyzing of seismic 
data. 

The procedure of localization consists in the following. Let us assume that we have an array of data 
{xi, X2, X3, xm, •••> %n} and take the initial sampling of the fixed length M. Then by passing through 
all array of values with the "work window" of the fixed length M we can calculated the time series of the 



noisy parameters {A,(T,T + Mr), X l (T + t } T + (M + l)r), A 4 (T + (N - M - 1)t,T + (N - l)r)} and 
{Aj(T, T + Mr), A;(T + r, T + (M + l)r), A. ( (T + (N - M - l)r, T + (N - l)r)} 



Obviously, it is inadmissible to use both very large intervals (M =1000 and more points) and very short intervals 
(M =50 points and less) for definition of local parameters \i(t), Ai(t). In the first case the physical sense of the 
localization procedure is lost. On the other hand, it is impossible to carry out any plausible correlation analysis with 



III. LOCAL NOISY PARAMETERS 
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small intervals because of gross errors. Therefore there is necessity for finding the optimum length of the initial local 
interval or quantity M. 

To determine the optimal minimal local sampling we have used the data corresponding to the calm state of the 
Earth before the technogenic explosion (an underground nuclear explosion) . The calculation procedure consists in the 
following. We have taken the interval of 40 points as the starting point and have calculated all the low-order noise 
parameters Ai (i — 1,2,3) and Aj (j — 1,2) by Eqs. © and jTUjl - Then the interval was consistently increased by 
unit time segment r and the relaxation parameters Ai , Aj were calculated every time at the increase of the interval. 
As a result of this procedure executed for the calm state of the Earth we have established that all parameters take 
"steady" numerical values at the interval approximately equal to 150 points and more (see Fig. 1 for more details). 
Namely, from Fig. 1 one can see, that the parameters Ai and Ai take minimal values in their absolute quantity at 
this length of interval, the amplitude of value fluctuations of the parameter A2 gets lower and level off. It is again 
the evidence of reduction of the noise influence. Fluctuations of the parameters A2 and A3 also decrease, and the 
parameters themselves take steady values starting with the sampling of length ~ 150 points. Applying this procedure 
to other data of the calm state of the Earth, we find the same behavior of noise parameters Ai, Ai and detect the 
minimal interval of 150 points again. So, we choose the interval of such length as being optimal for accumulation of 
local statistics llfl. 



IV. DEFINITION OF RELAXATION PARAMETERS Xi(t) AND A,-(t) FOR EARTHQUAKES AND 

TECHNOGENIC EXPLOSIONS DATA 

It is well known that modern seismic devices allow one to derive different quantitative and qualitative data about 
the seismic state of the Earth. In this work we analyse three various weak local earthquakes (EQ's) in Jordan 
(1998) [EQ(1), EQ(2), EQ(3)], one strong earthquake in Turkey (summer, 1999) [EQ(T)] and three local underground 
technogenic explosions (TE) [TE(1), TE(2), TE(3)] with the length of registration from 10000 to 25000 points. In 
case of strong EQ its seismogram contains 65000 registered points. All these experimental data were courteously 
given by the Laboratory of Geophysics and Seismology (Amman, Jordan) . All data correspond to transverse seismic 
displacements. The real temporal step of digitization r between the registered points of seismic activity has the 
following values, viz, r = 0.02s for the EQ(T), and r = 0.01s for all others cases. We defined such characteristics as 
the maximal amplitude of signal fluctuations before "event" a±, the maximal amplitude of signal fluctuations during 
"event" 0,2, EQ (TE) power a^/ai, time interval till EQ (TE) To, continuance of "event" T; and finally the total 
time of signal recording T to tai directly from seismic data. These quantities are presented in Table I. They give a clear 
notion about the duration and power of investigated phenomena. One can see from Table I, approximately 4500 — 5700 
points are accounted for the visible part a wavelet. This number of the recorded points allows one to execute analysis 
by means of the local parameters Ai (t) and Aj (t) . 

The procedure of calculation of time-dependence for Xi(t) and Aj(t) was based on the following operations. An 
interval with M ~ 150 points is taken, and noisy parameters Ai, Aj are calculated for it with the help of Eqs. JJjJ, l|10fl . 
Then the operation of " stepwise shift to the right" at the interval of the fixed length M is executed, and parameters 
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Table I: Some characteristics of Technogenic Explosions [TE(1), TE(2), TE(3)] and Earthquakes [EQ(1), EQ(2), EQ(3), EQ(T)] 
obtained from seismic data: a\ is the maximal amplitude of signal oscillations before the "event", a 2 is the maximal amplitude 
of signal oscillations at the time of the "event", T is the time from the beginning of the signal registration to the beginning 
the "event", Tj is the "event" duration, T to tai — N is the total time of signal recording. 





ai 




ai 




02/01 


T 


T 


Ttotal 


TE(1) 


1.03 • 10" 


-3 


8.24 • 10 


-3 


8 


4091 


4500 


12500 


TE(2) 


1.5 • 10" 


3 


18.5 • 10 


-3 


12.3(3) 


3538 


4500 


10000 


TE(3) 


1.18 • 10" 


-3 


56.5 • 10 


-3 


47.88 


4091 


4850 


15000 


EQ(1) 


0.74 ■ 10" 


-3 


7.1 • 10" 


3 


9.59 


5682 


4500 


15000 


EQ(2) 


0.94 • 10" 


-3 


6.1 • 10- 


3 


6.49 


12308 


5770 25000 


EQ(3) 


0.59 • 10" 


-3 


14.1 • 10 


-3 


23.9 


3182 


5700 


12500 


EQ(T) 


9 ■ 10~ 3 




20 




2.2(2) • 10 3 


« 13500 




65000 



are computed again. These actions are executed, while the initial sampling X(t) will not be finished. As a result we 
obtain the following dependencies { A 4 (T, T + Mr) , A; (T + r, T + (M + 1 )r) , . . . , Ai (T + ( TV - M - l)r, T + ( N - l)r) } 
and { Ai (T, T + Mr) , A, (T + r, T + (M + 1 )r) , . . . , A 4 (T + ( N - M - 1 ) r, T + ( N - 1 )t) } . If the character of the noise 
in the investigated data change, some signal will appear or disappear, and it will be directly reflected in the behavior 
of the relaxation characteristics. 

The results of the above described procedure for the case of EQ and TE data are shown in Figs. 2, 3 and Table 
II. However, in order to check up the optimized length of the local interval, we calculated local parameters Xi(t) and 
Aj(t) (i = 1,2,3), (j = 1,2) at the local sampling with the length M = 100,200,250,300,350 and 400 points. It 
turned out that a large number of various noises in the behavior of Ai (t) and Aj (t) are superimposed on the carrying 
trajectory at M = 100 points (parameters has a gross errors). The line-shapes of Aj(i) and Aj(t) practically cease 
to change at the sampling the length M — 150 and more, M = 200, 250, .... Once again it is evidence that the local 
interval with the length M — 150 points is optimal for the analysis of strong, weak EQ's and TE's. 

The detailed analysis of the received results allows one to reveal the following features. 

Weak EQ's and local underground TE's (for TE(3) and EQ(3), see Figs. 2, 3): 

1. All the parameters Ai(£) take only negative values (Aj(i) < 0), whereas the noisy parameters Aj(t) can be both 
positive and negative. 

2. Noisy parameter X\(t). The absolute magnitude |Ai| increases sharply in its amplitude by a factor approximately 
equal to 4 — 13.3 during EQ (various for different EQ's), and then it returns to its initial state, the restoration time 
T\ 1 and the duration of "event" T/ are approximately equal for weak EQ, i. e. T\ 1 w Tj. 

The absolute magnitude |Ai| also exhibits an abrupt rise ~ 2.2 — 3.5 times higher for TE's. However, it returns 
quickly to its normal level. The restoration time Ai (t) for TE is less than the duration of the " event" approximately 
by a factor of 2.5 — 3. 

3. Noisy parameter \2(t). The parameter A2(i) responds to the beginning of the "event" by an abrupt rise of its 
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value. It increases in amplitude both for TE's and for EQ's. The character of the noise changes during the "event". 
The parameter responds to the power of the " event" . 

4. Noisy parameter \^{t). This parameter always fluctuates near its numerical value — 1, and has an abrupt rise at 
the beginning the "event" in the form of separated spikes (see Figs. 3 and 4). The parameter keenly responds to the 
noise changes. 

5. Noisy parameter Ai(t). It fluctuates near zero before and after the "event" changing its sign at this time. The 
parameter increases sharply at the "event" and always (!) retains positive. Then it decays smoothly. Restoration 
time and the duration of the "event" Tj are approximately the same both for the EQ and TE. This parameter 
is very sensitive to the changes of the noise character. For example, the data analysis of the EQ(2) shows that the 
separate burst of the amplitude values of Ai(t) appears for ps 4000 points up to the beginning of the EQ, although 
such indicator was not visually observed in the initial seismic data [l7|. It may be the evidence of high prognostic 
property of this parameter for EQ's forecasting. 

6. Noisy parameter Az(t). Noise changes of the parameter A^(t) are observed during the "event" both during EQ's 
and TE's. From Figs. 2 and 3 we can see that Aa(t) has a distinctive negative depression during the "event". 

Strong EQ's (for EQ(T) see Fig. 4) Q: 

All parameters react keenly to the appearance of the signal in case of a strong EQ. Let us consider the behavior of 
the local parameters in this case in detail. 

7. Noisy parameter Ai (t) . The parameter demonstrates the presence of noise and takes negative values before the 
"event" [Region / in Fig. 4 f)]. As the "event" approaches, the amplitude of fluctuations decreases, and oscillations 
turn into negligible fluctuations near zero value for II, III and IV Regions (for more detail, see Fig. 4). 

8. Noisy parameter \2(t). A noise is also observed in the behavior of this parameter before the "event", and this 
parameter takes only negative values. However, its values decrease sharply in absolute magnitude and begin to take 
values near zero with the appearance of an EQ signal (Fig. 4 b). 

9. Noisy parameter As(i). This parameter takes only negative values at all times t. The negligible noise appears 
before the "event'. As the "event" approaches the noise increases in amplitude by the factor 3 — 5. The amplitude of 
oscillations begins to decrease in Region IV (see Fig. 4 c). 

10. Noisy parameters Ax(t) and Aa(t). Parameters fluctuate, taking positive values mainly before the "event" 
(Region I). With the beginning of "event" values of parameter increase greatly in absolute magnitude approaching 
zero. At the beginning of Region II both parameters change their sign from positive to negative. Then for Regions 
II, III and IV we see only right line with A\(t), A2(i) — > on the scales of Figs. 4 d) and e). However, the parameter 
A2(t) has faintly visible fluctuations for Region IV characterizing the final EQ phase [see Figs. 4 d) and e)]. 

So, all parameters are very sensitive to the approach of a strong EQ. A sharp change in their behavior is appreciable 
before strong fluctuations of the Earth's surface far off 10000 points (~ 3.5 minutes). 
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Table II: Characteristics of local parameters Xi(t) and Ai(t): T\ 1 is the parameter relaxation time, Ta 1 is the relaxation 

time of Ai(t), t\ is the relaxation time for exponential attenuation of X\(t), ta is the relaxation time for exponential attenuation 
of Ai(t). 





Ao (units of r 1 ) 


A A (units of r" 1 ) 


T A 


ta 


A (units of r 2 ) 


AA (units of r~ 2 ) 


T A 




AA/Ao 


AA/Ao 


Ti/T x 


TE(1) 


-0.13 


-0.425 


1800 


90 


0.02 


0.28 


4500 


45 


3.27 


14 


2.5 


TE(2) 


-0.15 


-0.34 


1440 


100 


0.02 


0.28 


4032 


55 


2.26(6) 


14 


3.125 


TE(3) 


-0.15 


-0.33 


1870 


80 


0.002 


0.42 


4850 


45 


2.2 


210 


2.59 


EQ(1) 


-0.045 


-0.6 


4500 


170 


0.005 


0.43 


4500 


110 


13.3(3) 


86 


1 


EQ(2) 


-0.1 


-0.5 


5770 


170 


0.001 


0.35 


5770 


130 


5 


350 


1 


EQ(3) 


-0.12 


-0.47 


5700 


210 


0.01 


0.35 


5700 


180 


3.91(6) 


35 


1 



V. SIMPLE EXPONENTIAL MODEL FOR TIME LOCAL BEHAVIOR OF NOISY RELAXATION 

PARAMETERS X x {t) AND Ai(t) 

As can be seen from Figs. 2 — 4 all relaxation parameters are very sensitive to the beginning of EQ and TE. The 
behavior of parameters X\(t) and Ai(i) in weak local EQ's and local TE's is of great interest. These parameters are 
initial in our calculations [see Eqs. @, i|l(J|) ]. The analysis has shown that parameters oscillate near average values Ao 
and Ao, correspondingly, before and after the oscillations visible on seismograms. The results of the behavior of these 
parameters for TE(3) and EQ(3) are presented in Figs. 2 a)-d) and Fig. 3 a)- d). However, a sudden rise by factors 
AAo and AAo is always observed in the behavior of these parameters at the enhancement or the appearance of the 
signal (it is seen at the beginning of the EQ or the TE in seismogram data). Furthermore, the continuous attenuation 
occurs. Over all this time these parameters have a well-defined pronounced trend. Such behavior of Ai(t) and Ai(t) 
give us a possibility of modelling the time dependence of these parameters by some simple mathematical functions. 
The fitting procedure showed that the time behavior of these parameters can be well approximated by the following 
simple time dependencies: 

Ai (t) = A + AA ■ exp f-*zJ?£ \ • H(t - T ), (46) 



Ai(f) = A + AA ■ exp {-^3} ' - r o), (47) 

where H(t) is the Heaviside function, T\ and Ta are the relaxation times of Xi(t) and Ai(i), correspondingly. The 
time To is the same in Eqs. (I46II and l|47|l for parameters Xi(t) and Ai(t) (see Table HJ. The numerical values of the 
variables, included in Eqs. (|46|l . I|47J) were defined for EQ's and TE's by comparison of localization results with these 
equations (see Fig. 5). Numerical values of parameters are presented in Table ITTl 

So, it proved that the restoration of these parameters to their steady values occurs according to the exponential law. 
As can be seen from Fig. 5 this description best suits for Ai(t) and Ai(i) of weak EQ's. The foregoing estimations 
strengthen fully our resume 2 and 5 of Section llVI 
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The results presented in the last three columns of Table II might be of interest for readers. As the quantities 
AAo/Ao and AAo/Ao show the rise of value of the corresponding parameter at TE and EQ. Finally, the ratio between 
the "event" duration Ti and the relaxation time T\ discovers a remarkable distinction between TE's and weak EQ's. 

VI. CONCLUSION 

In this work universal method for investigating non-stationary, unsteady and non-Markov random processes in 
discrete systems is suggested. This universality is achieved by combining the opportunities of microscopic, mesoscopic 
and macroscopic descriptions for complex systems. This method allows one to find and to analyze fast, slow and 
super-slow processes. To investigate super-slow processes we propose to use the model of a " quasi-Brownian particle" . 
The wavelet-transformation of the initial time series can be used for this purpose. This method helps to analyze and 
differentiate similar signals of different origin. Theoretical investigations have been realized by means of two methods 
supplementing each other: the statistical theory of discrete non-Markov stochastic processes 2] and the local noisy 
parameters. The application of the last method gives a possibility to study non-stationary and unsteady processes 
with alternation and superimposition of different modes. The correlation between the time scales characteristic of 
different modes may be different. However, at accurate realization the localization procedure allows one to separate 
the noise and the signals and to carry out their quantitative and informative analysis. 

Another important advantage of this method is the possibility to operate it in "real time" regime, i.e. it can be 
put into practice immediately at getting the data, that is of great practical value. 

The developed approach has been tested for strong and weak EQ's data and nuclear underground TE's. As a result 
we have obtained of the following results. 

The time-behavior of the local relaxation parameters can be described by simple model relaxation functions. The 
temporal relaxation of parameters Xi(t) and Ai(i) in weak EQ's and TE's after the beginning of the "event" occurs 
according to the exponential law. However, the restoration and duration of the events are practically the same in 
case of weak EQ's. The restoration time of parameters Ai(i) and Ai(t) in the case of TE's differs noticeably from the 
duration of the event. So, this approach can be useful in recognition these two different seismic phenomena. From 
the analysis of strong EQ data one can see that the behavior of all parameters changes greatly long before the EQ. 
For example, such change for the EQ(T) presented here occurs ~ 3.5 minutes before the main event. This change 
of Aj(t) and Aj(i) obtained for strong EQ's opens the possibility for a more accurate registration of the beginning of 
changes of the parameters before the visual wavelet and real EQ's. 

We are sure that the suggested method can be very useful for the study of a wide class of random discrete processes in 
real complex systems of live and of lifeless things: in cardiology, physiology, neurophysiology, biophysics of membranes 
and seismology. 
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Appendix A: THREE FORMS OF THE QUASI-OPERATOR I 

Equations of motion of a random variable x with the use of the derivative of the three different forms are 
represented here. 

Using evolution operator, we can write the equation of motion for a discrete case as following 

dx(t)/dt = i£x(t). (Al) 

However, there is a possibility of application of the time derivative d/dt — ► A/ At in three different forms: 
1. "Right" derivative (with decurrent difference in numerator) 

Ax{t) _ x(t + r) -x(t) _ 1 
At ~ t ~~ T 

with Liouville's quasioperator of the following form: 



U(t + r,t)x(t) (A2) 



C(t,T) = --[U(t + T,t)-l}; (A3) 

T 

2. "Left" derivative (with ascending difference in numerator) 

A*W = *{t)-*(t-r) = X (t)-U-Ht,t-r) X (t) = 1 ^ t 

At T T T 

with the Liouvillian of the next form: 

C(t,r) = --[l-U- 1 (t,t-T)}; (A5) 
r 

3. "Central" derivative (with central difference in numerator) 

Ax(t) _ x(t + t) - x(t - t) _ x(t + t) - x(t) x(t) - x(t - t) 
At _ 2r ~ 2r 2r 

= ^[^t + r.iJ-P-^M-TM). (A6) 
Then the quasioperator £ takes the following form: 

£&t) = -7T P(t + r, t) - U- 1 (t, t-r)}. (A7) 
It 

In the calculations and the analysis presented in this work we have used the derivative of the first form [see Eqs. 

1, Si. 
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Figure 1: The definition of the optimal size for a local "working window". Calculations were made for the data of a steady 
state of the Earth. As a result we find that size N ~ 150 points is more optimal. The relative stability of all parameters is 
observed. Inset f) shows seismic signals for the Earth's steady state (before underground TE). 
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Figure 2: The calculated time behavior of noisy relaxation parameters and A, for the technogenic explosion TE(3) [figures 
a)-e)], the signal of which is presented at inset f). 
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Figure 3: The time behavior of the noise relaxation parameters A; and A; for the weak EQ(3) [see figures a)-e)], the signal of 
which is also presented at inset f). 
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Figure 4: The time behavior of the noisy relaxation parameters A; and Aj [figures a)-e)] for the strong EQ(T), the signal of 
which is presented in figure f). The following regions: I - a calm state of the Earth's core, II - the state before the earthquake, 
III - the state during the earthquake and, finally, IV - the state after the event are divided by the vertical dotted lines. 




Figure 5: The time behavior of the first two noisy relaxation parameters Ai, Ai for TE(3) [figures a) and b)] and for the weak 
EQ(1) [figures c) and d)]. Solid lines show fitted functions 1 Kit and 14711 with corresponding parameters, presented in Table II. 
One can note the typical exponential restoration of the parameters with the beginning of the event. 



